typeset using JPSJ.sty <ver.l.Ob> 



0^ 
u 

< 

(N 



Charge Kondo effect toward a non-Fermi-liquid fixed point 
in the orbitally degenerate exchange model 
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We show that a Kondo-type model with an orbital degeneracy has a new non-Fermi-liquid 
fixed point. Near the fixed point the spin degrees of freedom are completely quenched, and the 
residual charge degrees of freedom lead to the multi-channel Kondo effect. Anomalous behavior 
appears in electric and thermal properties, but the magnetic susceptibility should show the local 
Fermi-liquid behavior. The non- Fermi-liquid fixed point becomes unstable against perturbations 
breaking the particle-hole symmetry. We derive these results using the third-order scaling for a 
spherically symmetric model with a fictitious spin. In contrast to the Coqblin-Schrieffer model, 
the present model respects different time-reversal properties of multipole operators. 

KEYWORDS: orbital degeneracy, time-reversal property, multipole, multi-channel Kondo effect, non-Fermi liquid, 
CeBe, LaMnOa 
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Orbital dynamics has recently attracted much interest 
because of new experimental progress in heavy-fermion 
systems such as CeBg^"^^ and transition metal oxides 
such as LaMnOa."*^^' In such systems, the localized 
ground multiplet has an orbital degeneracy n in addi- 
tion to the Kramers degeneracy. Thus orbital tensors 
as well as spin ones constitute observables.^"^-* The sim- 
plest model to describe the combined spin and orbital 
degrees of freedom is the SU(2n) symmetric spin-orbit 
model. In this model, however, the different time- 
reversal properties of the spin and orbital tensors are 
hidden in the excessively high symmetry. 

One can naturally ask how the Kondo effect occurs 
in the presence of the orbital degeneracy, and how the 
anisotropy in the combined spin-orbit space influences 
the physical property. In this context we refer to the pi- 
oneer work of Hirst who extended the Coqblin-Schrieffer 
(CS) model for the realistic atomic level structure. 
Pavarini and Andreani recently studied the stability of 
the local Fermi liquid in the CS model against the sym- 
metry breaking interaction found by Hirst, and con- 
cluded that it is stable. In order to include the realistic 
crystalline-electric-field (CEF) effect, the present authors 
considered more general exchange interaction between / 
and conduction electrons based on the point-group anal- 
ysis, and discussed how the time-reversal property of op- 
erators is reflected in the Kondo effect. In general, 
operators describing the spin and/or the orbital degrees 
of freedom are classified according to their definite time- 
reversal properties. The operators which are even under 
the time reversal are called electric, and those which are 
odd are called magnetic. The distinction between mag- 
netic and electric characters should be important in in- 
terpreting and predicting experimental results including 
dilution effects. 

In this Letter we propose a new mechanism to realize 
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the non-Fermi-liquid fixed point by combination of spin 
and orbital degrees of freedom. Near the fixed point 
only the spin degrees of freedom are compensated, and 
the residual charge degrees of freedom leads to the multi- 
channel Kondo effect. This non- Fermi-liquid fixed point 
requires the presence of particle-hole symmetry, and the 
anomalous behavior appears only in electric properties. 
Namely the magnetic susceptibility behaves as if the sys- 
tem is a local Fermi liquid. 

Let us start with our choice of the basis set and ex- 
plain their time-reversal property. The ground multiplet 
is written as \la) = 1, 2, ■ • ■ n, cr =t, j), where a repre- 
sents the Kramers pair with arbitrary strong spin-orbit 
interaction. The time-reversal operation 6 is defined as 



\la) = {-lf'^-'^\l-<j) 



(1) 



The conventional Pauli matrices tr" acting only on the 
spin state are odd under time reversal: 



(a = x,y,z). 



(2) 



We define the orbital operators which act only on 
the orbital state I. Their time-reversal properties are 
represented as 



1), 



(3) 



due to the antilinearity of the time-reversal operation. 
Thus, in the Hermitian representation the pure imagi- 
nary orbital operators are odd under time reversal, and 
have the magnetic character. In the simplest case of 
n = 2, the operator is magnetic and and are 
electric. In the case of n = 3, three of eight orbital oper- 
ators are magnetic and five are electric. In the absence 
of the spin-orbit interaction, magnetic orbital operators 
describe the orbital moment of electrons. 

In order to facilitate the scaling analysis, we map the 
states onto |m) where m represents the z-component of 
a fictitious spin j ^ n — 1/2 (half integer). An example 
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of mapping for n = 3 is given by 



|1T) 

|3i) 



+ 5/2), |2i) 
-1/2), |2T) 



+ 3/2), |3T) 
-3/2), lU) 



+ 1/2), 
-5/2)<4) 



We assume that the Hamiltonian has the spherical sym- 
metry in the new basis set. Then wc introduce the spher- 
ical tensor operators^^' which have a simple time- 
reversal property: 

^-1 = (_l)p[TW]t, (5) 

where each operator represents 2Pth magnetic (electric) 
miiltipole for odd (even) p with the component q. This 
property under time reversal merely reflects the odd 
character of the first rank operator, the angular momen- 
tum. 

The set of the spherical tensors with p = 1, 2, • • • , 2j 
together with the unit matrix completely describes the 
2n X 2n kinds of transition between local electron states. 
Then any spherically symmetric interaction among op- 
erators between / and conduction electrons is described 
in terms of the spherical tensors as 



Hint =^2 J, E ctT(rtc|/)[T(f)]t(/|, 

p=l q=-p 

-ip) 



(6) 



where Jp is the coupling constant, Tq is a matrix repre- 
sentation of the tensor operator with dimension 2n and 
summation over quantum numbers of each state is im- 
plied by the matrix multiplication. Here the spherical 
tensors are orthonormalized as 



Tr(TW[Tj'']t) = VV- 



(7) 



In general, the angular momentum j of conduction elec- 
trons is not necessarily the same as the localized elec- 
tron. However, we assume the same j for simplic- 
ity. In the case of Jp = Jcs for all p, one recovers the 
Coqblin-Schrieffer model with the SU(2n) symmetry. It 
should be emphasized that the spherical symmetry SU(2) 
in eq.(6) is adopted just for technical convenience. Ob- 
viously the actual system with the CEF is not invariant 
under the infinitesimal rotation. The simplified model 
presented above is sufficient for our purpose of analyzing 
the consequence of different time-reversal properties be- 
tween electric and magnetic degrees of freedom. Because 
the time-reversal property of operators is conserved un- 
der the discrete rotation of the point-group, the exten- 
sion of the model is straightforward. 

A general formula is available for constructing the scal- 
ing equation in terms of the structure constants of the 
algebra for the related matrices. ^^'^^^ The commutation 
rule is given by 



(8) 



where the Greek a, (3 and 7 are abbreviation of (p, g), 
{p',q') and {p",q"), respectively. We adopt the com- 
pletely antisymmetric structure constant which is given 



by 



p p p 



(9) 

where the parenthesis including three columns is the 
Wigner's 3j symbol. The symmetric function S with 
respect to p's is defined as 



Sg]p„ = V(2p+l)(2y + l)(2p" + l) 



p p' p" 
j j j 



(10) 

with p + p' + p" =odd. The bracket including three 
columns represent the 6j symbol. 

Once the commutation rule is given, we can easily ob- 
tain the scaling equations up to third order: 



d 



^ 7 _ g{2) I ^(3) 

'^p — yp ^ yp 



(11) 



2p + 
2J, 



-5^J,.V(5jJ^„)', (12) 



4'^ = 2^E^.^'E(^V)^ (13) 

pi pit 

where the coupling constants are measured in units of 
the bare cut-off energy Dq, and ^ = \5i{D/Dq). 

Let us simplify the analysis by putting Jp = J„i for odd 
p and Jp = Je for even p. Namely we pay attention only 
to the difference in time-reversal property of operators. 
Then eqs. (11)-(13) are reduced to 

^ = -2(1 - J„)(A« + A« J2), (14) 

— Je = 2A^J^ J(, [(J^ + Jg) — 2J„ 
where ^'s are defined as 

odd 



2p + 



^E(^i^V)'. (p = odd). 



p'p' 

even 



4^^ = ^E(^^;V')^ (P=odd), 



(15) 



(16) 



(17) 



odd even 



^^'^ = ^EE(4^V')'. (P = even). (18) 
p' p" 

The summation over p' and p" can be performed ex- 
plicitly and j4's are independent of p. For half-integer 
(integer) j, we obtain 

^W = l(2j + 3), Q(2i-1)), (19) 
4^) = i(2i-l), f7(2j + 3)V (20) 



^0) = 







(j = 1/2), 



(21) 



l;(2j + l) (otherwise). 

One can confirm that the scaling equation for the SU(2n) 
OS model is reproduced with J^ = J^ = Jgg. 

From eqs. (14) and (15), one can identify the following 
two fixed points for all n > 2: 
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(i) (ii) (j,;,j*)=(i,o). 

Clearly, the fixed point (i) has the SU(2n) symmetry and 
its local Fermi-liquid property has been studied well.^^-' 
The finite value of the coupling constant at the fixed 
point is an artifact of the third-order scaling and should 
be replaced by infinity in the exact scaling. 

The trajectory of the scaling equations is shown in Fig. 
1 for the case of j = 5/2 (n = 3). Note that Jm is driven 
to strong coupling no matter how small is Jg. On the 
other hand, Jg can retain the zero initial value even if 
Jm is finite. This aspect is related to the particle-hole 
symmetry in the system as will be discussed shortly. 
J 

»(i) 




Fig. 1. The trajectory of the scahng equations for j = 5/2 (3- 
fold orbital degeneracy). There are two fixed points {J^,J*)- 
(i) (1,1), and (ii) (1,0). 



For the SU(27i) CS model with 2n = 2j + 
1, the characteristic energy is given by = 

J^g^^"'' exp[— l/(2nJcs)]- In the absence of Je, on the 
other hand, the alternative energy scale appears: = 



l/(n+l) 



exp[— l/(n-f 1) J,„]. Figure 2 shows the solution 



of the third-order scaling equation for j = 5/2 {n — 3). 
Here the initial coupling constant is fixed at Jmo = 0.05 
and the ratio a = Je/ Jm is varied. It is emphasized 
that Je remains small over a substantial energy range 
although Jm is already in the strong-coupling regime. 
Hence one can identify two different characteristic ener- 
gies in the case of a <C 1. The two characteristic energies 
merge into a single Kondo energy as a approaches unity. 

Now we turn to discussing the property of the new 
fixed point (ii) with J* = 0. If one performs the particle- 



hole transformation 



'km 



the cou- 



pling constants Jp for the transformed system are related 
to the original ones as 



J n 



(22) 



which results from the anticommutation of fermion op- 
erators and the parity of spherical tensors. Thus the 
particle-hole symmetry is realized only for Jp — Q with 
p even. The SU(2n) CS model with n > 1 breaks this 
symmetry since the model necessarily has Je ^ 0. 

Let us discuss how the charge Kondo effect arises, as- 
suming that the initial coupling constants are already 
in the strong-coupling regime. The localized spin and 
one-particle or hole excitation from the Fermi sea form 




Fig. 2. The renormalization evolution for j = 5/2. The initial 
coupling constants are fixed at J^o = 0.05 with various values 
of the ratio a = J^o / JmO ■ 



the singlet bound state. We represent this state as (-)-, 0) 
where -|- means the binding with a particle excitation and 
the singlet. There is the particle-hole conjugate state 
(—,0) where a hole is bound to the localized spin. The 
energies of the two states are different by the amount Jg. 
If one reduces the magnitude of Je to zero with Jm fixed 
in the strong-coupling regime, the ground state with the 
two-fold degeneracy in the charge degrees of freedom, 
(±,0), is realized. Therefore, at the fixed point (ii) the 
spin degrees of freedom is completely compensated, but 
the charge degrees of freedom (isospin) still remain. 

We now construct the effective interaction around the 
fixed point (ii). In the lowest order of t^/Jm, with t 
being the hopping of innermost conduction electrons, 
there appear interaction terms which connect the two- 
fold degenerate ground states. These terms accompany 



-k'm 



, or equivalently (—1)-' ™'c\^c] 



and the 



matrix element (+, 0|cj.^cj,,_^| — , 0) is finite. Since the 
singlet state is isotropic, the pair creation with any of 
m has the same magnitude in the interaction. Thus the 
effective interaction around the fixed point (ii) is given 

by 



(23) 



where the "isospin" operators for conduction electrons 
are defined as 



J 



fc 

\ / , ^km^k' -m^ L^mJ 



kk' 

and for localized electrons as 

-•4 



| + ,0)(+,0|-|-,0)(-,0| 



X+ = | + ,0)(-,0|, I- = [I+] 



(24) 
(25) 

(26) 
(27) 



The effective interaction is equivalent to that of the 
multi-channel Kondo model with the overscreening. The 
multiple channels come from the orbital degeneracy. 
Here we have included the weak particle-hole asymmetry 
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as the second term with h in eq.(23), which is of the order 

of Jg. Just like the multi-channel Kondo model which is 
unstable in the presence of magnetic field, the electric 
coupling Je is the relevant perturbation in the present 
model. 

Since the ground state has a degeneracy with J,, — 0, 
the system exhibits an anomalous behavior with the 
fractional entropy as in the conventional multi-channel 
Kondo effect. ■^'^^ The effective Hamiltonian is meaning- 
ful only in the case where the renormalized coupling 
becomes strong while Je still is in the weak-coupling 
regime. In such a case only the electric response Xp{'^) = 
(Ti^^(a;)Ti^)) with even p shows non-Fermi-liquid be- 
havior. On the contrary Xpi'-^) with odd p shows a 
Fermi-liquid behavior because the spin degrees of free- 
dom is quenched. We are now exploring quantitatively 
the unique behavior of the system by using the Wil- 
son's numerical renormalization-group (NRG) method, 
and some anomalous properties, such as the fractional 
entropy, have been confirmed. ^^-^ 

In summary, we have investigated a multipolar ex- 
change model with the rotational symmetry. Respect- 
ing the different time- reversal properties of operators, 
we find a new fixed point where a non-Fermi liquid ap- 
pears as a result of the Kondo effect in the charge sec- 
tor. Although the fixed point is unstable against any 
perturbation breaking the particle-hole symmetry, the 
non-Fermi-liquid behavior may be observed over a sub- 
stantial energy range in the case of Je <C Jm- Thus, more 
quantitative investigation of the model is highly desired. 
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